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Atmospheric refraction is responsible for the bending of light-rays in the atmosphere. It is a result 
of the continuous decrease in the refractive index of the air as a function of altitude. A well-known 
consequence of this phenomenon is the apparently elliptic shape of the setting or rising Sun (or 
full- Moon). In the present paper we systematically investigate this phenomenon in a standard 
atmosphere. Theoretical and numerical calculations are compared with experimental data. The 
asymmetric rim of the Sun is computed as a function of its inclination angle, observational height 
and meteorological conditions characterized by pressure, temperature and lapse-rate. We reveal and 
illustrate some extreme and highly unusual situations. 



I. ATMOSPHERIC REFRACTIONS 

Atmospheric refraction is between the first scientifi- 
cally approached phenomenon. This effect is responsible 
for the apparent scintillation of stars, mirages, the puz- 
zling and spectacular green flash, difference between ap- 
parent and real position of stars or the asymmetric rim of 
the Sun near the horizon. The refractive index of dry air 
is very close to 1. Its small dependence on temperature 
and pressure leads to a refractive index gradient in the 
atmosphere. Although this gradient is very small, due to 
the large distances travelled by the light-rays in the at- 
mosphere, it can result observable bending or dispersion. 

Atmospheric refraction is usually divided into three 
categories: astronomical, terrestrial and geodesic. Ter- 
restrial refraction appears when both the object and ob- 
server are within the Earth's atmosphere. This refraction 
is responsible for ordinary mirages, and has been exten- 
sively studied Q|. The first one reporting and account- 
ing for mirages was Aristotle in Meteorologica 0. Cor- 
rectly, he concluded that dense air layers can act as mir- 
rors, and considered this effect responsible for mirages. 
(A complete historic and bibliographic study of mirages 
can be found on the splendid web-page of A.T.Young 
P|.) Geodesic refraction is a special case of terrestrial 
refraction where both the object and observer are at low 
altitudes. A well-known example for this is surveying. 
We speak about astronomical refraction when a terres- 
trial observer detects ray-bending effects for the light- 
rays coming from objects outside the Earth's atmosphere. 
This refraction is responsible for the difference between 
the real and apparent position of stars near the horizon, 
the green flash or the asymmetric rim of the setting (or 
rising) Sun. In the present paper we discuss in detail 
the flatness of the Sun or full-Moon near the horizon. 
Particularly, we are interested in the extreme values this 
flatness can take for standard atmospheric profiles, and 
the dependence of the flatness on the altitude of the ob- 
server and meteorological conditions (temperature and 
pressure). A similar study, considering fixed atmospheric 



conditions was recently published by Thomas and Joseph 
Q. Here we plan a more complete analysis, considering 
three different theoretical methods, experiments, com- 
puter simulations and illustration of the results with a 
collection of pictures and movies. 

The spectacular nature of sunset or sunrise makes the 
phenomenon attractive to students. After our experience 
it can effectively be used to exemplify refraction phenom- 
ena and exercising the principles of geometrical optics. 
Our paper is structured as follows. First, we introduce 
the atmosphere model, and the refractive index profile 
is given. Then, we present three different methods for 
computing the ray-path in this atmosphere. Results com- 
puted for different meteorological conditions and observa- 
tion altitude will close the theoretical part. Experimen- 
tally, we analyze and compare with theory the measured 
flatness for some sunsets photographed or video-filmed 
by us. Finally, we discuss some extreme and unusual 
conditions which are exemplified by pictures and movies 
on a web-page supporting this study. From this web- 
page one can also freely download a computer-program 
written by us, which visualizes the sunset for arbitrary 
meteorological conditions and observation altitude. 

We emphasize here, that our study is restricted to stan- 
dard atmospheric conditions with smooth temperature 
and pressure profile. Non-standard, but often encoun- 
tered atmospheric profiles leads to non-standard distor- 
tions of the solar rim, and are not considered within this 
study. 

II. THE OPTICAL ATMOSPHERE MODEL 
(REFRACTIVE INDEX AS A FUNCTION OF 
ALTITUDE) 

The atmosphere of the Earth is composed mainly of N% 
(79%) and 2 (20%), and extends up to a few hundreds 
kilometer height. The refractive index of air is very close 
to 1, and depends slightly on its pressure and temper- 
ature following Edlen's semi-empirical law M, valid for 
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dry air: 



1 + 10" 6 (776.2 + 4.36 x 10"V) 



T ' 
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In the above formula v is the wave-number of the light in 
cm , P is the pressure in kPa and T the temperature 
in K. Since the pressure and temperature varies within 
the atmosphere, we get a refractive index gradient which 
is responsible for atmospheric refractions. Although the 
variations in n are quite small, the large distances trav- 
elled by light-rays in the atmosphere makes refraction 
effects observable and sometimes important. 





Ionosphere 

100- 




1 1 

p p f 

Pessure (hPa) 


80- 

Mesosphere a 

— 60- 
*r 

3 




Stratosphere 40 " 




— IU 


Tropopause 

20 




- 100 


Troposphere 
^^^^^^^^^^^ 





Earth 



Temperature (°C) 

FIG. 1. Distinct layers of the atmosphere, and the standard 
temperature and pressure profile. 



is the U.S. Standard Atmosphere, established in 1953 
and re-actualized in 1976 @. The U.S. Standard At- 
mosphere consists of single profiles, representing the ide- 
alized steady-state atmosphere for moderate solar activ- 
ity. The listed parameters include temperature, pressure, 
density, gravitational acceleration, mean particle speed, 
mean collision frequency, mean free path, etc. as a func- 
tion of altitude. In our study we consider an atmosphere 
model with a spherical symmetry, all relevant physical 
quantities (temperature and pressure) varying only as a 
function of altitude. We then calculate the refraction in- 
dex of air as a function of altitude as follows: 
(i) The temperature profile in the troposphere is linearly 
decreasing with a A = 6.5K/km lapse-rate, as suggested 
by the U.S. Standard Atmosphere. Following again the 
U.S. Standard Atmosphere we consider the temperature 
constant within the tropopause. Although the temper- 
ature increases as a function of altitude in the strato- 
sphere, due to the rarefied air (small pressure) the re- 
fractive index is close to 1, and for the sake of simplicity 
we assume the temperature as constant in this region, 
too. (We checked that this approximation is fully jus- 
tified.) Above z s = 50km height we assume that the 
refractive index is 1, and do not calculate it anymore by 
Edlen's formula. Up to the top of the stratosphere the 
temperature profile is therefore presumed as 



T(z) 
T(z) = To - 



= T - Xz for z < z t 
Xz t for z t < z < z s , 
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As a function of altitude several layers with different 
physical properties are distinguishable. The lowest layer 
extending from sea-level to approximative z t — 14fcm 
height is called the troposphere. This is the region where 
the wheatear takes place, i.e. the region of rising and 
falling packets of air. In this layer the air pressure drops 
drastically, at the top of the troposphere being only 10% 
of the value measured at sea- level. In the troposphere the 
temperature decreases almost linearly with the altitude 
(Fig. 1). A thin buffer zone between the troposphere and 
the next layer (the stratosphere) is called the tropopause. 
Within the tropopause the temperature is in good ap- 
proximation constant. The stratosphere extends from the 
altitude of approximately 18km up to z s = 50km (Fig. 1). 
Within the stratosphere the air flow is mostly horizontal. 
In the upper part of the stratosphere we have the ozone 
layer. Inside the stratosphere the pressure decreases fur- 
ther with the altitude, but surprisingly the temperature 
increases with height. Above the stratosphere we find 
the mesosphere and the ionosphere. In these regions the 
air is very rare and the temperature profile is shown in 
Fig. 1. 

From the viewpoint of atmospheric refractions only the 
first two layers, the troposphere and the stratosphere are 
important. In the upper layers of the atmosphere the air 
is so rarefied, that the refractive index can be considered 
1 within a good approximation. 

An accepted and widely used model for the atmosphere 



where z is the altitude from sea-level, and To the tem- 
perature at sea-level. 

(ii) For the pressure profile we use a barometric formula 
in which we take into account the variation of tempera- 
ture with altitude. Considering a vertical slice of air with 
thickness dz, the variation of pressure within this slice is 
due to the hydrostatic pressure 



dP = —p(z) g{z) dz, 



(4) 



where g(z) is the gravitational acceleration and p(z) the 
density of air at height z: 



*> - v 



Nm P(z)m 



V 



kT(z) 



(5) 



(We denoted by m the mass of one molecule, N the num- 
ber of molecules in volume V, T{z) the temperature at 
height z and k the Boltzmann constant.) Since we fo- 
cus on the troposphere and the stratosphere only, the 
z height is small in comparison with the radius of the 
Earth (Re ~ 6378fcm), thus the g gravitational acceler- 
ation can be considered constant. We can write thus 



dP = - 



P(z)mgdz 
kT(z) ' 



(6) 



which is a separable differential equation for P(z). Inte- 
grating this between a height zq where the pressure is Pq 
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and an arbitrary height z, we get the desired barometric 
formula: 



P(z) = P exp (-^ jf 



dz' 



(7) 



Considering z Q = (sea-level) and using the temperature 
profile given by equation (||) and (||) a simple calculus 
leads us to 



P(z) = Po 
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for z < z t , and 
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mg(z - z t ) 
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3) in Edlen's 
stratosphere 



for z t < z < z s . Plugging this and (t 
formula (|l]), we get for the troposphere and 
the variation of the refractive index as a function of al- 
titude. As stated before, for altitudes higher than z s we 
can merely use n = 1. The refractive index profile as 
a function of altitude calculated in this manner for the 
parameters of the U.S. Standard Atmosphere is plotted 
in Fig. 2. 
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FIG. 2. The used refractive index profile for To 
P = latm, A = 6.5K/km and v = 20000cm" 1 . 
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III. TRAJECTORY OF A LIGHT-RAY IN THE 
ATMOSPHERE 

We present now three different methods for computing 
the path of a light-ray in our optical atmosphere model. 
We proved that the results given by these methods are 
the same, justifying our forthcoming theoretical consider- 
ations. None of these methods is purely analytical, they 
all make use of numerical calculations to derive the light- 
ray trajectory. In the following we briefly describe these 
methods and sketch how one can compute the deviation 
angle due to atmospheric refraction for light-rays coming 
from distant sources. In order to avoid the phenomenon 
of dispersion let us first consider monochromatic light. 



A. The Integral method 

This method closely follows the one discussed by Smart 
0. To illustrate the method we use the geometry from 
Fig. 3. 




FIG. 3. Real (5V) and apparent (S) position of a distant 
source S, as observed from the point M within the Earth's 
inner atmosphere. 

An observer in M (at altitude Zq) detects the light- 
source placed in S r (altitude z'). Let the arc S r M be the 
presumed path of the light-ray between the S r source and 
the observer in M. The observer detects the light ray in 
the MS direction, which is tangent to the ray path in 
M. The source will be positioned by the observer in S 
(SO and S r O have radial directions, thus the light-rays in 
these directions would not bend within our optical atmo- 
sphere model). We denote by 5 the apparent inclination 
angle, characterizing the direction of the S image. Let 
Sdev be the deviation angle of MS relative to the real 
MS r direction of the source. The meaning of the 9 and 
r angles are obvious from the figure. We are now inter- 
ested to compute Sdev as a function of 6. Presuming that 
z' >> zo (i.e. the source is very far from the Earth), 
we can approximate MS, MS r and OS by z' . Some el- 
ementary geometry will convince us that the following 
approximations are justified: 
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In order to get the desired (W«(<5) dependence we need 
9 r as a function of <5. 

Let us follow now a light-ray approaching the Earth 
and let us consider the atmosphere stratified in infinites- 
imally thin layers of thickness Ar, with slightly different 
refractive indices (Fig. 4). 
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FIG. 4. Trajectory of a light-ray in a layered optical atmo- 
sphere model. 

A useful relation between the initial incident angle ao 
(in layer with refractive index n ) and a later incident 
angle a k (for a layer with refractive index rife) can be 
derived. Using Snell's law and the notations from Fig. 4, 
we can write: 



smao ni 



sin a\ 
A0i = ai 



"o 
i 



This leads to: 



sinai = — sinao (cos A9i + sin A9i cota'i), 
m 



(12) 
(13) 
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Similarly 
sin a 2 



ni 

"2 



sinai(cos A0 2 + sin A# 2 cota 2 ), (15) 



and combining with (|l~4|) it leads to: 

sin a 2 = — sin ao(cos A9\ + sin A9\ cot a\ ) x 

(cos A6» 2 + sin A9 2 cot a' 2 ). (16) 
Generalizing the above equation we get: 

sinafc = — sinao (cos A9\ + sin A0\ cot a!) x ... x 

rife 

(cos A6>fe + sin A6» fc cot a fe ). (17) 

Since the A0 4 angles are small, the following approxima- 
tions are justified: 
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(18) 
(19) 

(20) 



Using these approximations and the fact that Ar << n 
we get: 
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Presuming now r /r k w 1, we get: 

rono sin ao 



smafe = 



r-fcrife 



Following the geometry from Fig. 4 we get 

Ar Ar 1 
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r,; r, 
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(21) 



(22) 



(23) 



which in the Ar — > and at — a(r) continuous limit 
yields: 



d9 = 



dr 
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r I i i 

sin a(r) 



(24) 



Denoting by am the final incident angle at the observer 
and using (|22l) we get: 



dS = 



dr 



1 



(25) 



ri(r) 



n(r M ) 



- 1 



sin^ a M r M 

By using the fact that sin 2 a a/ = cos 2 (5 we can finally 
determine the angle 9 r as a function of <5 by integrating 
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The above integral can be only numerically computed. 
There is of course a singularity at z = Zq, which can be 
eliminated by a Gauss-Chebisev expansion, or by adjust- 
ing the step in the numerical integration in the vicinity 
of the singularity (the method followed by us) . Comput- 
ing numerically 9 r as a function of 5 and by getting from 
( |l0| ) 9 as a function of 6, we are able now to calculate 
from ([ll]) the desired Sdev(5) dependence. 

It is important to mention here that applying the 
method for 6 < is not straightforward. One must first 
find in this case the closest point C of the light-ray tra- 
jectory relative to the Earth's surface. This can be done 
by using the fact that the trajectory is symmetric in the 
vicinity of this point. Then, we decompose the trajec- 
tory in two parts, the first part is from M to C, and 
the second from C to S r . The light-ray can be followed 
now by the presented methods both on the CS r and CM 
segments. The deviation angle can be also computed. 
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B. Using the Fermat principle 

The Fermat principle states that light travels between 
two point along that path which requires the least time, 
as compared to other nearby paths. We can of course re- 
formulate the Fermat principle by using the optical path 
instead of time. If the trajectory of a light-ray travelling 
in the X — O — Y plane is described by the y = y(x) curve 
(see the geometry in Fig. 5), we have that: 

rPf in {x t ,yf) 

s = n[x, y(x)]ds = 

J Pini(xo,ya) 



n[x, y{x)] a/1 + y'(x) 2 dx, 



(27) 



should have a local minima. In the above formula 
n[x,y(x)] denotes the refractive index of the medium in 
point with coordinates a;, y(x). In our case n has a spher- 
ical symmetry, depending on the z = \/ x 2 + y{x) 2 — Re 
altitude only. 
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FIG. 5. Geometry and notations for the method based on 
the Fermat principle. 



One can write then ( |27| ) as: 

f[x,y(x),y'(x))dx, 



(28) 



with: 



f[x, y(x), y'{x)] = n[^Jx 2 + y(x) 2 - R E ] y/l + y'(x) 2 . 

(29) 

We are looking now for the y{x) function that minimizes 
s. If we fix the two points Pini and Pfin between which 
the light travels, the minima of s leads to a classical vari- 
ational problem 



5s = 0, 



(30) 



with: 



S[y(x)] U=x = o, 
S[y(x)] \ x = Xf = 0. 

The solution of this problem is well-known 
by the second-order differential equation: 
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It is straightforward to show that: 
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By simple algebra we obtain from this a second-order 
differential equation for the y(x) equation describing the 
trajectory of the light-ray: 

/// \ dn [1 + y'(x) 2 } 

y [x) = — [x,y(x)\ = x 

dz n[x,y(x)]y' x 2 + y(x) 2 

{y(x)[l + y'(x) 2 } - y'{x)[y(x)y'(x) + x]}. (36) 

Taking the n{z) refractive index profile from our opti- 
cal atmosphere model, equation ( pq ) can be numerically 
integrated. We start from a Pini{xo,yo) point and con- 
sider a tan S — y'(xo) initial derivative. The angle 5 
will be the apparent inclination angle of the light-ray 
in Pini (the point where the observer is presumed). The 
y(x) trajectory can be computed by numerically integrat- 
ing ( [36] ) with the y(xo) = yo and y'(xo) = tan S initial 
conditions. We construct thus the y(x) trajectory from 
point-to-point, up to an altitude z > z s , where n(z) = 1 
can be presumed. The derivative of y(x) at this point 
will determine the 5 r angle of the light-ray at the bor- 
der of the optical atmosphere. For altitudes higher than 
this, the trajectory of the light-ray is presumed recti- 
linear. The deviation angle is then simply approximated 
as Sdev — 5 — Sr, and the desired 5(5dev) dependence can 
be numerically computed. 



C. The simulation method 

As its name suggests this is simply a computer sim- 
ulation method in which we follow the light-ray by seg- 
ments of infinitesimally small and fixed-length lines. The 
angle between two elementary line is given by estimat- 
ing the refractive index at their end-points and by using 
Snell's law. The trajectory of the light-ray is then con- 
structed starting from the observer with an initial angle 
<5, and computing the ray-path from point-to-point until 
it leaves the optical atmosphere and the refractive index 
can be taken as 1. At this point the direction of the 
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light-ray determines the final angle S r , and the deviation 
is computed as Sdev = S — S r . Resembling the method 
based on the Fermat principle, this simulation also works 
smoothly for arbitrary positive or negative 8 values. Due 
to the fact that for an elementary step the changes in the 
refractive index is rather small, it is crucial to work with 
the maximal precision offered by the computing environ- 
ment. 

We mention here that all three methods give identi- 
cal results, proving their applicability and our theoretical 
considerations. As an example, for the parameters of the 
U.S. Standard Atmosphere, results for the Sdev deviation 
angle versus the apparent inclination angle, S, is plot- 
ted on Fig.6 (T = 10°C, P = latm = 101.3fcPa, v = 
20000cto~ 1 wave-number for the light and A = 6.5K/km 
lapse-rate) . 
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S(degrees) 

FIG. 6. Deviation angle Sdev between the apparent and real 
location of a point-like source far from Earth as a function of 
the 8 apparent inclination angle. (To = 10°C, Po = latm, 
standard A = 6.5K/km lapse-rate and observation from 
sea-level) 

These results are in excellent agreement with the one 
given by Thomas and Joseph [3 and the report of the 
U.S. Naval Observatory (1993) |§. From Fig. 6 we learn 
that the deviation angle is usually quite small, and be- 
comes important only when viewing objects in the vicin- 
ity of the optical horizon. The deviation angle increases 
sharply for very small inclination angles. This effect is re- 
sponsible for the flattened shape of the rising (or setting) 
Sun and Moon, and also for the fact that these heavenly 
objects appear more flattened at their bottom. It is also 
interesting to note here that the angular extent of the 
Sun and full-Moon is 32 arcmin, and the maximum de- 
viation obtained for standard conditions is 34.5 arcmin. 
This leads us to the observation that the Sun or Moon 
is visible even when in reality it is below the geometrical 
horizon. 



IV. COMPUTING THE FLATNESS 

Once we numerically determined the Sdev deviation an- 
gle as a function of the 5 apparent inclination angle, it 
is easy to characterize the apparently flat rim of the set- 
ting Sun. Atmospheric refraction influences only the A„ 
vertical angular extent of the Sun (or Moon), which be- 
comes thus smaller than the A^ horizontal extent. The 
asymmetry ratio (or the flatness) for the rim of the Sun 
can be described by the 



(37) 



ratio (of course a G [1, oo)). Since both the Sun and the 
full-Moon are normally visible under a Ao — 32 arcmin 
angular extent, we have A/j = Ao. The value of A„ can 
be derived after numerically computing the <5 apparent 
inclination angle as a function of the 5 r = 8 + Sdev real 
inclination, i.e. S = F(S r ). If the apparent inclination 
angle for the bottom of the Sun is So, corresponding to a 
S r o real inclination angle, than 



A„(<$o) = F(<y r o+Ao)-<y , 



and we get: 



a{S ) 



A 



F(S r0 + A ) - S 



(38) 



(39) 



For fixed observation altitude and meteorological condi- 
tions the a c maximal possible flatness corresponds to the 
situation when the bottom of the Sun touches the hori- 
zon. This happens for a So = 5 C critical inclination angle. 

By decreasing the S angle in small steps, and apply- 
ing the previously described methods for computing the 
light-ray trajectory, both the S — F(S r ) function and S c 
is numerically determined. In our calculations we have 
chosen to decrease S in steps of 0.01°. 



V. RESULTS 

Applying the methods presented in the previous sec- 
tions and our optical atmosphere model, we systemati- 
cally computed the a c asymmetry ratio for the rim of 
the setting Sun as a function of observational altitude 
and meteorological conditions parameterized by pressure 
and temperature. We have also shown that the value of 
a is insensitive to details of the considered optical atmo- 
sphere model, proving the stability of our results. In the 
following we detail our findings. For all the calculations 
we considered monochromatic light with wave-number of 
v = 20000cm" 1 , corresponding to the green color. If oth- 
erwise not specified we considered the height of the tro- 
posphere z t — 14km and the height of the stratosphere 
as z s = 50km. 
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A. Stability regarding the considered optical 
atmosphere model 

As discussed in section II., for the construction of 
the refractive index profile we borrowed results from the 
U.S. Standard Atmosphere model. The main parameters 
needed by us is the Zt height of the troposphere, the z a 
height of the stratosphere and the A lapse-rate. 

1.2267 i ■ 1 ■ 1 ■ 1 
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FIG. 7. Maximal observable flatness, a c as a function of 
the z t height of the troposphere. (To = 0°C, To = latm, 
normal lapse-rate and observation from sea-level) 
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FIG. 8. Maximal observable flatness, a c as a function of 
the A lapse-rate in the troposphere. (T = 0°C, Po = latm, 
and observation from sea-level) 

It is evident, that the exact value of z s does not much 
influence our results, since in the stratosphere the refrac- 
tive index is already very close to 1. As a first step we 
studied thus the influence of Zt on the a c asymmetry ra- 
tio. We considered normal conditions with To = 0°C, 
T = latm, observations at sea-level (z — 0) and 
A = 6.5K/km standard lapse-rate. As illustrated on 
Fig. 7, the value of z t (in a reasonable range) has no 
significant influence (note the scale on the vertical axis). 
The value of the lapse-rate has already a more noticeable 
effect on a c (Fig. 8), however this variation is also quite 
small for the practically important fluctuations around 
the standard A = 6.5K/km value. We conclude thus, 
that our results are quite stable regarding the details of 



the considered optical atmosphere model. 

B. Asymmetry as a function of the inclination angle 

We computed the observable a flatness as function 
of the apparent inclination angle of the Sun's bottom. 
Results for T = 0°C, T = latm = lOlMPa, A = 
Q.5K/km and observations at sea-level are presented on 
Fig. 9. As emphasized before, the asymmetric rim of 
the setting (or rising) Sun becomes evident only for very 
small Sq values, when the Sun is close to the horizon. For 
these normal parameters we get that a c is around 1.2. 

1.3 | . , . , . , . 1 
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FIG. 9. Observable flatness, a, as a function of the So incli- 
nation angle of the bottom of the Sun. (To = 0°C, To = latm, 
observation from sea-level and standard lapse-rate) 

C. Flatness as a function of observation height 

Let us presume now that the observer is at height z 
above the sea-level, and there is no obstacle in the direc- 
tion of the horizon, which is at sea-level (i.e. we are above 
a vast ocean) . It is obvious that from higher altitude the 
5 C critical angle will be smaller (becomes negative) and 
the deviation angle increases more sharply in the neigh- 
borhood of S c . This leads us immediately to the conjec- 
ture that the observed flatness should be larger. Taking 
the To = 0°C and To = latm normal atmospheric condi- 
tions at sea-level, the standard A = Q.5K/km, lapse-rate 
we can compute the a c maximal observable flatness as a 
function of observation height. Results in this sense are 
plotted in Fig. 10. As expected, a c increases with zq. 
For these normal conditions we get that from the top of 
a 5km height mountain we would observe an a c w 1.5 
flatness and from a commercial flight at 10km height at 
sunset we would detect an a c w 1.7 maximal flatness. 
For altitudes above 30km, or observations made from a 
space-shuttle we can get extreme values for a c up to 2.5. 
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FIG. 10. Maximal observable flatness, c 
the zo observation height. (To = 0°C, Po 
dard A = 6.5K/km lapse-rate) 
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as a function of FIG. 12. Maximal observable flatness, a c as a function of 

: latm and stan- the atmospheric pressure at sea- level. (To = 0°C, observation 
at sea-level and standard A = 6.5K/km lapse-rate) 



D. Influence of temperature 



VI. EXPERIMENTS 



The Tq temperature measured at sea-level influences in 
an important manner the observed flatness. Computing 
the temperature dependence of a c for observations at 
sea-level, normal Po = latm atmospheric pressure and 
standard lapse-rate we get the values presented in Fig. 11. 
While for a To = 30°C temperature a c is around 1.1 for 
T = — 40° C it becomes 1.3, and increases more and more 
sharply for lower temperatures (arctic conditions). 
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FIG. 11. Maximal observable flatness, a c as a function of 
the air temperature measured at sea-level. (Po = latm, ob- 
servations at sea-level and standard A = 6.5K/km lapse-rate) 



E. Influence of atmospheric pressure 

Increasing the P pressure at sea-level results the in- 
crease of the observable a c flatness. For the reasonable 
and measurable Po values, T = 0°C, observation at sea- 
level and standard lapse-rate the variation is almost lin- 
ear (Fig. 12). For an extreme Po = 125fcPa pressure one 
can detect an asymmetry ratio of 1.3. 



By simple experiments it is relatively easy to measure 
the flatness of the setting (or rising) Sun. We considered 
photo and video experiments, and analyzed the pictures 
as a function of the inclination angle of the Sun. With 
an appropriate filter and calibrated eye-piece, telescope 
observations were also possible. In order to obtain us- 
able pictures with a nicely visible rim we had to ensure 
a properly adjusted light intensity. Appropriate filters, 
specific atmospheric conditions and usually small 6q in- 
clination angles lead us to usable pictures. When the 
Sun is close to the horizon, the inclination angle of the 
Sun can be directly determined by analyzing the taken 
picture. We can use the fact that the horizontal angular 
extent of the Sun always corresponds to 32arcmin. The 
distance of the Sun from the horizon can then be com- 
pared with the horizontal extent of the setting Sun and 
the inclination angle results directly from the picture. 
Making experiments with the setting or rising Moon is 
more complicated, since one needs a full-Moon for this, a 
work during the night and we have to deal with pictures 
where the horizon is not clearly visible. 

First, we have taken several series of pictures in South- 
Bend (Indiana, altitude 100m) both in winter and late- 
spring, studying very different temperature conditions. 
As an immediate confirmation of our theoretical results, 
from these pictures it was obvious that for the same incli- 
nation angle the flatness is bigger in winter, i.e. for lower 
T values. Two series of digitized and appropriately en- 
hanced pictures are visible on the web-page dedicated to 
this study |10|. Results from these pictures, in compar- 
ison with the expected (computed) flatness is presented 
on Fig. 13. The first set of pictures is a sunrise in winter. 
The mean temperature during the sunrise was — 5°C, and 
the atmospheric pressure was 103fcPa. It is important to 
note that during this sunrise the temperature remained 
approximately constant. The second set of picture was 
made in late-spring for a sunset. During this sunset the 
mean-temperature was 12 C, the atmospheric pressure 
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98kPa and the temperature dropped detectably during 
the time the pictures were made. Results for the flatness 
as a function of the Sun's inclination angle (bottom of 
the rim) is plotted with filled circles and triangles for the 
winter and spring series, respectively. 
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FIG. 13. Flatness as a function of inclination angle. Com- 
parison between experimental results and theory for a sun- 
rise and a sunset photo sequence. Filled circles are results 
from sunrise pictures taken in winter with To = — 5°C and 
Pa = 103kPa. Filled triangles correspond for pictures taken 
in May with < T >= 12°C and P = 98kPa. Observations 
were made at zo = 100m height with the horizon roughly at 
the same altitude. Theoretical results for the corresponding 
atmospheric conditions and A = 6.5K/km standard lapse-rate 
is plotted by continuous and dashed lines for the winter and 
spring conditions, respectively. 

The theoretical curves were constructed for the men- 
tioned mean temperature, pressure, observation height of 
100m with the horizon at the same altitude and a stan- 
dard lapse-rate. These results are plotted on Fig. 13 by 
a continuous and dashed line for the winter and spring 
conditions, respectively. As observable from Fig. 13 for 
the winter conditions the measured a values are in ac- 
ceptable agreement with the one given by our theory. 
However, for the spring series the expected a values are 
higher, and decreasing slower as a function of the <5o in- 
clination angle, than the measured data. A reason to 
account for this sharper trend is the decreasing tempera- 
ture measured during the sunset. Taking this effect into 
account in the calculations, would definitely result in a 
trend closer to the observed one. This temperature vari- 
ation does not account however, for the constantly lower 
values measured for a. The only arguments we can give 
in this sense is that probably the atmosphere at the time 
of this measurement had a refractive index profile differ- 
ent from the one proposed in our model, behaving in a 
non-standard manner. 

A second set of experiments were realized by video- 
filming sunsets in Cluj (Romania). The obtained contin- 
uous set of picture allowed us to follow-up more precisely 
the flatness as a function of the inclination angle. The 
inclination angle was calculated from the images by the 
same method as in photographs, i.e. by comparing the 



height between the bottom of the Sun and horizon with 
the horizontal extent of the setting Sun. Since Cluj is 
not a flat region like South-Bend, we had to be care- 
ful in choosing the observation point, and to determine 
also the altitude of the optical horizon. Results from 
video-recording in comparison with theoretical expecta- 
tions (corresponding to the appropriate atmospheric con- 
ditions) are plotted in Fig. 14. For this measurement a 
quite fair agreement between theoretical and experimen- 
tal data is achieved. Since the lapse-rate for the theoret- 
ical prediction was taken from the standard atmosphere 
model, the slightly greater a values calculated by us are 
explainable. 
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FIG. 14. Experimental (dots) and theoretical results 
(dashed line) for the flatness as a function of the inclination 
angle from a sunset video-recording. (To = 7°C Po = 108fcPa, 
observation at zo = 400m, horizon at zo = 300m and a pre- 
sumed standard A = 6.5K/km lapse-rate) 



VII. PICTURES AND MOVIES FROM THE 
INTERNET 

Sunset or sunrise (moon-set or moon-rise) is usually a 
spectacular phenomena, and thus it is a favorite theme 
for professional and amateur photographers. The In- 
ternet is full of beautiful and useful pictures in this 
sense. Many of these pictures are taken in extreme con- 
ditions (arctic environment, airplanes, space-shuttle or 
high mountains), offering us an excellent possibility to 
check our calculations under these conditions, too. More- 
over, on the web there are also interesting, freely down- 
loadable or public domain movies exemplifying how the 
observed flatness increases in the neighborhood of the 
horizon. 

We performed an extensive search on the Internet and 
collected non-copyrighted materials about sunset, sun- 
rise, moon-set and moon-rise. Since their presentation 
in the context of this paper is impossible we classified 
and stored them on the web-page [[l0| dedicated to this 
study. The interested reader can browse this collection 
and convince himself or herself that these pictures sup- 
port our theoretical results for the estimated flatness. In 







agreement with our expectations we found that for a sun- 
rise (or sunset) viewed from a space-shuttle a c should be 
of the order 2 — 2.5, and in arctic environment a c in- 
creases up to values of 1.3. For most of the everyday, 
usual and low-altitude photos one finds a c » 1.1 — 1.2. 
Pictures taken from commercial airplane yield a c ss 1.5. 
We also found a picture-series for a sunset over the ocean, 
where the rim of the Sun is nicely visible. The scenario 
presented in this photo is again in agreement with our 
theoretical predictions. 

The movies stored on our web-page will also convince 
the reader about the sharp variation of a as a function 
of the Sq deviation angle. A public domain movie, show- 
ing a moon-set viewed from the space-shuttle illustrates 
observations from high altitudes. 

VIII. THE GREEN FLASH 

Although the explanation and the study of the green 
flash is not the aim of the present paper, we briefly dis- 
cuss here how this phenomenon can be understood and 
studied through our methods. 

From equation (Q) it is obvious that atmospheric re- 
fraction depends on the frequency (color) of the light- 
ray, a phenomenon called dispersion. Differently colored 
light-rays coming from a source emitting a continuous 
spectra will suffer different deviations. More strongly 
will bend the light-rays with bigger frequencies, leading 
to a separation of the colors in the observed Sun. Since 
the high-frequency visible components (corresponding to 
violet and blue colors) are strongly scattered by the at- 
mosphere, the green component reaching directly the 
observer will suffer the strongest bending. When the 
Sun disappears bellow the horizon this component will 
be observed lastly, leading to a green-flash on the hori- 
zon. Since for standard atmospheric conditions the dis- 
persion from atmospheric refraction is tinny, the effect 
is hardly observable for usual and everyday conditions. 
Green flashes become observable exactly under those con- 
ditions under which the flatness is more accentuated. 
This means high observation altitudes, low temperature 
and high pressure. As an illustration of this, we have just 
learned that for airplane pilots it is a quite usual phe- 
nomenon. Green flashes are often seen in non-standard 
optical atmosphere, when mirages appear. This is gen- 
erally the case when layers of air with strongly different 
temperatures are in contact. 

For a nice presentation and discussion on green flashes 
the interested reader should consult the excellent home- 
page of A.T. Young g. 

IX. SUNSET SIMULATION PROGRAM 

We also created a computer-code which simulates the 
sunset on the computer-screen. By using the discussed 



simulation method, the program computes and visual- 
izes the rim of the Sun during a sunset. After fixing 
the atmospheric conditions (temperature and pressure) 
and observation height in the menu, one can follow up 
how a sunset might look like in our optical atmosphere. 
The program runs under Windows environment and one 
can freely download the executable from the web-page 
accompanying this study |l(J ■ 

X. CONCLUSIONS 

Atmospheric refraction is responsible for the asymmet- 
ric rim of the setting (or rising) Sun. Three different 
methods yielding the same results were presented here to 
compute the path of a light-ray in an optical atmosphere 
model, where the refraction index varies continuously as 
a function of altitude. By determining the deviation an- 
gle between the apparent and real inclination of a point- 
like light-source which is at a large distance from the 
Earth, we were able to compute the a asymmetry ra- 
tio for the rim of the setting Sun. We investigated a as 
a function of the inclination angle, observation altitude 
and atmospheric conditions. We found that the maximal 
flatness obtained in the vicinity of the horizon increases 
as a function of observation altitude and pressure, and 
decreases as the temperature increases. We found that 
a is rather insensitive to the fine-details of the consid- 
ered optical atmosphere model, which makes our results 
robust. The maximal flatness observable under normal 
conditions at sea-level is around 1.2. At extremely low 
temperatures (— 40°C) one can observe values up to 1.3 
and for high altitude observations (space-shuttle) one can 
get an asymmetry ratio of 2.5. By simple experiments 
and pictures from the Internet we illustrate and prove 
our theoretical predictions. The methods presented here 
can be effectively used to study atmospheric refractions 
in non-standard atmospheric conditions as well. In this 
manner one can study mirages or the green flash phe- 
nomenon. A freeware computer program created by us 
and downloadable from the web-page accompanying this 
study simulates the sunset for arbitrary atmospheric con- 
ditions and observation altitude. 
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